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ABSTRACT 



The diffuse galactic gamma- ray spectrum meas ured 
by the EGRET experiment l|Hunter et al.lll997j) are 
interpreted within a scenario in which cosmic rays 
(CRs) are injected by three different kind of sources, 
(i) supernovae (SN) which explode into the interstel- 
lar medium (ISM), (ii) Red Supergiants (RSG), and 
(iii) Wolf-Rayet stars (WR) , w here the two latter ex- 
plode into their pre-SN winds (|Biermann et al.l20011 
iSina et al.ll2f)nl . 

Key words: cosmic ray sources; gamma-rays; diffuse; 
galactic . 



1. INTRODUCTION 



Conventional models of diffuse galactic gamma- 
ray production are based on three main processes. 
Gamma-rays are produced through the decay of tt" 
as secondary particles of hadronic collisions of CRs 
with the ISM, such as proton-hydrogen or proton- 
helium collisions, in brcmsstrahlung processes of CR 
electrons with the ISM, and through inverse Comp- 
ton scattering of CR electrons with interstellar radia- 
tion fields. These models can explain a wide range of 
observations like the energy spectra below ^ 1 GeV 
or the integrated flux from the outer parts of our 
Galaxy. Nonetheless the diffuse galactic gamma-rays 



Table 1. Interaction spectra for the different type of 
supernovae. ii^kneo and Ecut-oS are in GeV. 71 and 
72 are the spectral indices = ^qE^ ) below and 
above i^knoc- Z is the charge of the nucleus. 
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Figure 1. The log of the mean energy of the CR par- 
ticles versus the number of shock crossings for Vsh = 
0.01c. Starting from the bottom for 5° , 25° , 65° , 80° 
and 85° respectively. We see the difference in the 
energy gain of CR particles for the almost perpen- 
dicular case co mpared to sm aller shock inclinations. 
An effect that ^Jokivii\198^) pointed out as well. 



observed by the EGRET experiment from the inner 
Galaxy above w 1 GeV exceeds by about 60 % the 
intensity predicted by these calculations (the mea- 
sured spectrum is too hard). 

In this contribution we shall present first results on 
the diffuse gamma-ray production expected in the 
model byEicrmanii et al. (20^. In this model, in 
addition to the aforementioned processes, interac- 
tions of hadrons of the relatively hard CR injection 
spectrum (see Tab. 01 with the pre-SN winds are con- 
sidered as further possible sources of diffuse galactic 
gamma-rays (Bicrmann ct al. 2001). 

In section 2 a Monte Carlo code is developed to 
model particle acceleration in non-relativistic near 
parallel and highly oblique shock configurations, in- 
cluding cross-field diffusion, with application to RSG 
and WR winds. In section 3 we review the adopted 
model of diffuse galactic gamma ray production and 
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introduce all necessary input parameters. In sec- 
tion 4 and 5, we present our results, compared to 
EGRET, CASA-MIA and KASCADE data. In sec- 
tion 6, we summarize our findings and discuss some 
further improvements in our models. 



2. SOURCES OF COSMIC RAYS 



iJokiniil l)1987|) investigated the rate of the energy gain 
and the maximum energy in non-relativistic shocks, 
which can be attained in given conditions, such as the 
effect of a highly oblique magnetic field to the scat- 
tering of the CR particles. Briefly, he showed that if 
the perpendicular diffusion {k±) is much smaller than 
the parallel one CR particles can gain consider- 
able energy in quasi-perpendicular shocks compared 
to quasi-parallel ones. 

We have constructed a Monte Carlo code in order 
to simulate non-relativistic near parallel and oblique 
shocks with 0° < t/i < 90° where tp is the angle be- 
tween the shock normal and the magnetic field, seen 
in the shock frame. The velocities of the upstream 
plasma flow are kept between 0.001c and 0.01c, which 
correspond to astrophysical environments with non- 
relativistic shocks such as in WR winds, RSG winds, 
etc. In order to investigate the behavior of the parti- 
cle scattering, between the quasi-parallel and oblique 
shock configurations, a cross-field diffusion is also al- 
lowed. Briefly, our Monte Carlo simulations in com- 
parison to the Jokipii's theoretical work show the 
same trend (see Fig. Furthermore, for realis- 
tic environments, such as RSG and WR stars, our 
calculations suggest that WR winds could provide 
large enough radii for protons to reach enough en- 
ergy in the available time and space, and we claim 
that the maximum energy (for Z=l) reached could 
be - lO^"^ eV, in about 10^° sec (w lO^yrs), which 
corresponds to lO^^cm (« lOpc). On the other hand 
in RSG winds we may find the conditions (smaller 
radii and corresponding time) where the maximum 
energy could be about 3 x 10^"* eV. There is further 
work under way which calculates the real time versus 
energy, attainable in highly oblique non-relativistic 
shocks found in WR and RSG stars and investigates 
in detail the constraints put in these models, fol- 



lowing 'l imits' that may apply (e.g. n 
l|MeH fc Bicrmann 20041 
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3. MODELS FOR GAMMA-RAY PRODUCTION 



[b] The gamma-ray spectrum is computed using 
a three-dimensional spatial distribution of matter 
and radiation fields in our Galaxy. In this model 
the ISM is composed of atomic (« 1 x IO^Mq), 
molecular (« 3 x lO^M©) and ionized hydrogen (w 
W^Mq) as well as 10 % He. Its spati a l distribution is 
taken from 'Laun hardt et al.l l)2002t) ; iKalberla et al.l 
l|1998r i. A logarithmic spiral arm model based on 
COBE/DIR BE data is als o incorporated (see Fig.|21) 
((Drimmel fc SDerge]ll200ll . The interstellar photon 
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Figure 2. Model of the ISM density based 
on COBE/DIRBE data Wrimmel & SveraelWoOl.) . 
The logarithmic spiral structure is superimposed on 
a purely radial distribut ion of molecular and ion- 
ized h ydrogen ifLaunhardt et al\\200^ : \Kalberla et al\ 
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Figure 3. Gamma-ray production spectra, see 
Eq. ^\). The spectra are shown in units of photons 
per atom, or crrr' in the case of inverse Compton. 



densities of interest have been approximated by three 
diluted blackbody spectra from the interstellar ra- 
diation field (ISRF) plus the 2 .7K microwave back- 
ground (CMB) (BloemeDill985H . 
CR spectra and their spatial distributions are an- 
other input for the calculations. We assume CR pro- 
tons and electrons to be radially distributed with a 
radial exponential scale a — b kpc. The assumed 
exponential scale heights hn ji are 0.5 kpc for pro - 
tons and 2 kpc for electrons ( Bloemenlll985l FigSGj) . 
The spectral indices are approximated by constants 
throughout the Galaxy and the normalization of the 
flux is taken from the review of IWiebel-Sooth et al.l 
EMI- 



The number of gamma-rays from a certain direction 
{l,h) with an energy per unit of time can then be 



Figure 4- Here the solution of lL{l,b) (Eq. 0) is 
shown for different combinations of CR and target 
particle types k,T. 



= J dL{l, b) J dEk 



calculated using 

dn~^{E-y, I, b) 
dtdE-ydn 

<Pk{Ek,r, I, b)ai{Ek,E^)nTir, I, b) 

where / dL{l, b) is the line-of-sight integral along the 
direction Ek and ^k{Ek,r,l,b) are energy and 
flux at position {r,l,b) of the CR particles of type 
k, ai{Ek, Ej) is the cross section of the gamma-ray 
production process i, and nxir, l,h) is the density of 
target particles of type T (either components of the 
ISM, photons of a radiation fields or particles of a 
SN-wind). In the case of IC it is also necessary to 
integrate over the energy distribution of the photon 
gas e. If we use the assumed independence of CR 
energy spectra from the position in our galaxy 



^k{Ek,r,l,h) 



$f(£;,)-<i>^(r,;,6), 



the integral can get factorized into one integral over 
energy and one over the line-of-sight: 



dnj{E^, I, b) 



dtdE^dfl 



Ie{E^) ■ lL{b,l), 



with 



and 



Ie{E^) 



J dEk-i>^iEk)a,iEk,E^ 



lL{b,l)^ dLnT{rJ,b)^l(r,l,b). 



(1) 



(2) 



The cross sections ai for Bremsstrahlung and inverse 
Compton effect have been numerically computed 
using the formulae given by J31umenthal & Gould 
l|l97Cll . For nuclear rea ctions they have b een sim- 
ulated with DPMJET3 ijEoesler et al.ll2n01t) . 
The independent solution for Ie and II are shown 
in Fig. El and Fig.H 



SPECTRAL ANALYSIS OF THE 
GAMMA-RAY EMISSION 



The emission from the inner Galaxy and from the 
outer Galaxy, predicted by our model calculations. 



energy (GeV) 
speL-lmm - -2 < b < 2 and 50 < I < 200 



B 



m ol all LOnlributions - 
pp- gamma- ISM - 
pH^-gamma-lSM 
pp-gamma-RSG - 
pHe - gamni a-RSG 




1000 10000 1000(H) le-(-06 

energy (GeV) 



Figure 5. Gamma-ray energy spectra ( complete 
model) for the directions of the galactic center (A) 
and anti-center (B). 



are shown in Fig.|ni(A) and (B), respectively, where 
the predicted spec tra are compared to EGRET 
l|Himter et alJl997^■ Whipij le (LeBohec et al. 2000), 
CASA-M IA ("Borion e et all mm and KASCADE 
(|Schatz e t al. 2003) data. The contributions to the 
gamma-ray spectrum from the three kind of super- 
novae, from Bremsstrahlung and from inverse Comp- 
ton scattering, as well as their sum, are explicitly 
indicated. The contribution to the diffuse gamma 
emission due to the collisions of CR protons with 
protons and helium in the interstellar medium ex- 
plains the spectrum above 0.1 GeV and below 1 GeV. 
Around 1 GeV the observed flux of diffuse galactic 
gamma-rays stems mainly from interactions of CRs 
injected by both ISM and RSG supernovae. Above 
1 GeV it is produced essentially by interactions of 
protons injected by exploding RSG with the protons 
in the local strong and enriched winds. Finally we 
remark that a change in the spectral index of the 
gamma-ray emission is evident from Fig. [S] (A) at 
photon energy above 1 GeV. Within our model the 
change in the gamma-ray spectral index arises nat- 
urally from the assumption that the two different 
regimes in the gamma-ray emission, corresponding to 
gamma-rays emitted by different supernovae types. 
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Figure 6. These plots only show the contribution to 
the gamma-ray flux predicted by conventional mod- 
els. The SN-wind component is not included. Mean 
flux for —2 < b < 2. Different contributions 
to the longitudinal distributions: Inverse Compton, 
Bremsstrahluny, Nuclear, Diffuse background, Sum. 
The EGRET data is shown as black dots. 



5. DIRECTIONAL ANALYSIS OF THE 
INTEGRATED FLUX 



The predicted longitude intensity distributions of dif- 
fuse galactic gamma radiation are shown in Fig. 
As can be seen in Fig. El (a) the used 3D models for 
ISM, photons and CR particles could be further op- 
timized to reproduce the EGRET data even better 
in this energy range. More importantly, it is appar- 
ent (Fig. El (b)-(d)) that for higher energies the flux 
coming from the Galactic center is getting more and 
more underestimated while the flux from the anti- 
center stays relatively well reproduced. The discrep- 
ancy between the calculation and the EGRET data 
supports the picture of an additional contribution 
from CR sources (SN) that are much more abundant 
in the inner Galaxy than in the outer. 



before they explode as supernovae. This wind 
material provides most of the grammage crossed by 
CR particles seen at Earth. We expect the data in 
the TeV range from MAGIC and HESS experiments 
to allow further insight or put constraints on the 
present model. More detailed studies are planned 
to calculate the variou s predictions of the model by 
iBiermann et al.l l)2001() and improve the 3D of the 
interstellar matter, cosmic ray, and radiation field 
distributions. 
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6. DISCUSSION AND OUTLOOK 



The presented model is very promising in explaining 
the measured diffuse gamma-ray flux. Also it seems 
to be capable to rep roduce the C/B ratio as well as 
the antiproton flux l|Sina et alJl2001|) . It is quite a 
plausible scenario in which the CR interact mostly 
in the environment close to their sources. Some 
stars, like RSG or WR, provide sufficient material, 
ejected as powerful winds at the end of their lives. 



